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Abstract 

Electrical pacing is a common procedure that is used in both experimental and clinical settings for 
studying and/or annihilating anatomical reentry. In a recent study [Comtois and Vinet, Chaos 12, 
903 (2002)], new ways to terminate the one-dimensional reentry using a simple protocol consisting 
of only two stimulations were discovered. The probability of annihilating the reentrant activity is 
much more probable by these new scenarios than by the usual local unidirectional block. This paper 
is an extension of the previous study in which the sensitivity of the new scenarios of annihilation 
to the pathway length is studied. It follows that reentry can be stopped over a limited interval of 
the pathway length and that increasing the length beyond the upper limit of this interval yields to 
a transition to sustained double-wave reentry. A similar dynamical mechanism, labeled alternans 
amplification, is found to be responsible for both behaviors. 

PACS numbers: 87.19.Hh, 05.45.-a 
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I. INTRODUCTION 



The picture of a fixed waveform traveling at constant speed around a ring of excitable 



tissue, is still a common representation of 
larly in reference to common atrial flutter 



unctionaljreentry in the clinical setting, particu- 
. However, the findings that complex 



reentries are possible even in a simple homogeneous one-dimensional ionic loop model and 
that their occurrence is dependent on the steepness of the restitution curve of the action 



in 



potential duration has altered the current understanding of the phenomena 
which reentry was postulated to remain stable and periodic as long as there was a minimal 
excitable gap ahead of the wavefront. These findings have also altered the thinking about 



the effect of antiarrhythmic drugs. 
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Overdrive pacing using a transvenously inserted catheter in the right atrium is a standard 
clinical procedure to interrupt atrial flutter. It is very successful, particularly when it is 
applied in conjunction with the administration of class I or III antiarrhythmic drug 3, 3 
The use of rapid pacing is likely to increase with the implant of permanent single or dual site 
stimulator for the prevention of atrial tachycardias Q]. However, the mechanism by which 
overdrive pacing interrupts reentry and the electrophysiological parameters of the reentry 
circuit that may determined an optimal choice of parameters for the pacing protocol are not 
understood. 

As a first step to improve the pacing algorithm, we have previously studied of a simple 
protoco 
reentry 

unidirectional block 



of stimulation consisting of two electrical stimuli applied in the pathway of a periodic 
I5I ]. New scenarios of reentry annihilation were identified, different from the classical 
18^ . which is still considered to be the most important mode of 
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termination. 

These alternative scenarios of reentry annihilation follow from a spatiotemporal process 
that we have called alternans amplification The first objective of this paper is to 

understand the effect of length of the reentry pathway on these scenarios of annihilation. 
We also show that beyond a critical length of the reentry pathway, alternans amplification 
induces a transition double-wave reentry instead of annihilation. 
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II. MODELS AND METHODS 



Results obtained with two different models are presented. The first model (ionic loop: 
IL) is a one-dimensional reaction diffusion system, using a cardiac ionic model to represent 
the transmembrane currents. The second model is an integral-delay equation (ID) based on 
the local properties of propagation and repolarization. 

A. Ionic loop model 

The well-known monodomain cable equation for an ID homogeneous excitable cardiac 
tissue embedded in an unbounded external medium of negligible resistivity is: 



where V is the transmembrane potential (mV), C m is the membrane capacitance (1 /iFcm~ 2 ), 
5* the surface-to- volume ratio (0.4 /im" 1 , assuming cylindrical cells with radius of 5 jum) and 
p is the intracellular resistivity (200 f2cm). I ion is given by a modified version of the Beeler- 
Reuter model (MBR) of the cardiac cell membrane, whose details and space-clamp dynamics 
are given in 

For each time step (At = 2//s), the system becomes a second order ordinary differential 
equation that is computed with a Galerkin finite element method projected on a linear basis 



function and a regular spatial mesh (Ax = 50 /zm)pj. The resulting tridiagonal linear 
system of equations is solved with a simplified LU decomposition method. The choice of 
At and Ax is motivated by the fact that depolarization is the stiffest part of the process. 
Programs were written in Fortran77 and ran on SGI workstations (Silicon Graphics). 

Following the generation of an action potential, two quantities are measured at each site 
to analyze the propagation: the activation time (T act (x)) and repolarization time (T rep (x)). 
T ac t{x) corresponds to the onset of the action potential and is defined as the moment at which 
dV/dt reaches its maximum during the upstroke of the action potential. T rep (x) is meant 
to indicate the time from which a new action potential can be generated by an incoming 
activation front or an external stimulus. A large set of simulations of sustained reentry in 
loops of different lengths has showed that an active propagating response was generated if 




(1) 
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stimulation was applied at least 30 ms after the —50 mV post-upstroke downcrossing in 
repolarization. Accordingly, this instant (i.e. 30 ms after —50 mV crossing) is taken as 
T rep (x). The action potential (A (x)) is defined as the time interval from T act (x) to T rep (x), 
such that T rep (x) = T act (x) +A(x). The diastolic interval (D (x)) associated to an activation 
occurring at T act (x) is defined as the time between the previous T rep (x) and current T act (x). 
With these definitions, a site is excitable if D (x) > 0. 



B. Integral-delay model 

The integral-delay model used in this study is an extension of a previous formulation that 
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was developed to describe sustained unidirectional propagation on the loop 0, 
A first relation A SP (D) gives the duration of the action potential as a function of D in 
the space-clamped configuration. If the nodes were disconnected from their neighbors, the 
repolarization time following an activation occurring at T act (x) would be: 

T r % (x) = T act (x) + A SP (D (a;)) 

provided that D > 0. If D < 0, the node is unexcitable. The actual repolarization time 
of a node at position x is expressed as a weighted average of T r s ep (x) over a symmetric 
neighborhood of length 2a, i.e.: 

/a 
w (u) T^Jp (x + u) du 
a 

where , w (x) = exp (— (3x 2 ) with = 0.5, (3 = 800 cm -2 , a = 0.5 cm), is the weighting 
function representing the effect of resistive coupling on the repolarization phase. (3 controls 
the spatial decay of the weighting function and d is a normalization coefficient such that 
J" w (u) du = 1. The calculation of T rep associated to one excitation is performed at 
each node at the next instance when it is stimulated by an incoming front or an external 
stimulus. At this moment, the Tf ep associated to the previous excitation of each point of 
the neighborhood are collected and averaged to produce T rep . In this way, a front whose 
propagation stops at some location still produces a continuous distribution of T rep around 
the region of block since T rep is an weighted average of the T^ e P of the sites excited by 
the front that is blocked and of those that were not reached by that front and still have 
the T^ ep (x) associated to their previous excitation. It provides at once a representation of 



the acceleration of repolarization of the excited cells induced by the load of those that are 
not excited, and of the prolongation of repolarization in those that were not excited by the 
electrotonic depolarization induces by the proximal excited cells. 

Once the T rep associated to the last excitation that we label T^p(x), has been computed, 
the diastolic interval associated to the current stimulation, which take place at the time 
Tact(x), is calculated as: 



D{x)=T act {x)-T l ; e ;\x) (2) 

If D > 0,the stimulus produces an action potential, which propagates on both side with a 
conduction time c t (D(x)), and reached the neighbouring nodes at the time T act (x ± Ax) = 
Tact (x) + Ax c t (D(x)). If D < 0, the point is not activated, and its T^ e P (x) is not changed. 

The integral-delay model was originally developed to represent the propagation of a 
single activation front during reentry, without external stimulation. In this context, T^^x) 
is always the repolarisation associated to the previous passage of the activation front and it 
can be written as T rep (x — L) = f " Q w (u) T^ e P (x — L + u) du = f_ a w (u) (T ac t (x — L + u) + 
A SP (D (x — L + u)) du. Similarly, T act (x) = T ac t (x — L) + c t (D(u))du. With these 
relations, eq. becomes: 

D (x) = Tact (x) - T rep (x - L) (3) 

pot 

= Tact (x — L)— w (u) Tact (x - L + u) du 



-a 

pL pa 

/ c t (D(u))du- / w(u) A SP (D(x - L + u))du 

JO J -a 



If the conduction time is constant J"" w (u) T ac t (x — L + u) du = T ac t (x — L) , and eq. El is 
equivalent to the version of the ID model introduced in jj. In fact this version neglects the 
effect of the delay of propagation in calculating the effect of coupling on repolarization. If 
w (u) is taken as a 5 function, which is equivalent to ignoring the effect of coupling, eq. El 
corresponds to the version of the integral-delay model of Courtemanche and al. 
The simulation of the ID model were performed using 

c t (D) = c tMn + 17.645e~ D / 2 - 79 + 12.40e" D / 15 - 34 (4) 
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where Cf > m m = 14.09 ms/cm, D is in ms and c< is in ms/cm, and 

A SP (D) = 42.23 + [Ao - + 30 ) 4 ' 92 (5) 

v > F 4 - 92 + (£> + 30) 4 - 92 



with A = 225.22 ms, B = 105.74 ms, r = 94.71 ms, and F = 59.38 ms jg [15j. These 
functions were obtained by fitting the data gathered from different regimes of propagation 
(free running periodic and QP reentry obtained with the IL model) . Computations were 
performed with a spatial discretization of Aa; = 50 /im, as in the ionic model. 



C. Protocol of stimulation 



For the IL model, I s um consists of a 2.5 ms current pulse applied over an interval of 450 
fim with amplitude of 60 /iA/cm 2 . Dual stimulations were applied on periodic reentry. The 
timing of stimuli was controlled by two parameters: Apsi, the time interval between the 
last activation at the center of the stimulated area and the onset of the first stimulus Si, 
Asi_s2 , the time interval between the onset of S\ and that of the second stimulus 5*2 

For the ID model, the value of D, defines the response of the nodes that are stimulated. 
If D < 0, the stimulus is applied in the refractory period and doest not have an effect. 
If D > 0, the stimulus depolarizes the tissue, defining A SP (D), and induces bi-directional 
propagation. As for the IL model the stimulation covers 450 /xm. 



III. RESULTS 



The post-stimuli dynamics are constrained by the steady states of the system. For loops 
longer than a minimum length L m [ n , sustained reentries are stable attractor of the system. 
These sustained reentries can be either periodic (period-1) or quasiperiodic (QP), and hold a 
single (SW), two (DW) or more traveling activation fronts. Table lists the stable solutions 
of both the ID and IL models for L < 3L m i n . The number and nature of the sustained 
reentries change with L and condition the outcomes of the stimulations. 

QP reentries are characterized by a spatial oscillation of D, with a wavelength A that 
is an irrational fraction of L (fig. ^\ and C). For both SW and DW reentries, there is an 
interval of L in which two different QP solutions coexist. These solutions, labeled mode-0 
and mode-1, have a similar structure but A mode_0 pa 3A mode_1 . In SW reentries, the passage of 
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Interval 


Reentry type 


L (cm) 


SW (single wave) 


DW(double wave) 


< L < L min = 12.8 






Lxmn ^ L < L^f^ e ^ = 15.2 

111111 II11I1 


QP, mode-0 




LZn ' 1 < L < L crit = 19.5 


QP,mode-0 
QP, mode-1 


— 


L C rit — L <C 2L m [ n = 25.6 


period- 1 




2L min < L < 2L^f n dc - 1 = 30.4 

— 111111 


period- 1 


QP, mode-0 


2L-f n dc " 1 < L < 2L CTit = 39 


period- 1 


QP, mode-0 
QP, mode-1 


2L cr n < i < 3L m j n 


period- 1 


period- 1 



TABLE I: Stable reentries of the IL model as a function of the length of the loop. 



each activation front is associated with a profile of A and D holding one (mode-0, panel A) 
or multiple (mode-1, panel C) maxima and minima over two turns. For DW QP reentries, 
L is twice the value for SW QP reentries, such that mode-0 has one maximum and one 
minimum over one turn. (fig. d^). Successive activations at each site alternate between long 
and short A and D values, except at a number of nodes corresponding to the boundaries 
from which the phase of the alternation is inverted. As illustrated by the time-course of 
T act (x) and T rep (a;)(fig. HJ right panels), the quasi-periodic nature of the propagation makes 
the position of the extrema and of the nodes to drift slowly in the direction inverse to the 
jropagation of the activation fronts. QP reentry is thus constituted by discordant alternans 
22.013 WJ th boundaries between short and long APD moving around the loop. 



Period-1 SW reentry can be annihilated by an isolated stimulus applied in the narrow 



vulnerable window in which the stimulation produces on 
to the well-known mechanism of unidirectional block 



V a retrograde front, corresponding 
12, (3 I3- m a previous paper, 



we have also described other modes of annihi 
that were induced by two successive stimuli 



ation as well as different transient dynamics 



3- 



These new modes of annihilation were 



compatible with experimental observations, and relevant to antiarrhythmic pacing therapy 
261 ] . However, that study was restricted to a specific length of the loop (L crit < L = 25 cm 



< 2L min ). In the following, we present a systematic investigation of the outcomes of double 
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FIG. 1: IL model: the mode-0 (panels A-B) and mode-1 DW reentry (panels C-D) coexisting at 
L = 34 cm. In this case, two fronts F\ and F2 travels along the loop. Left panels: the ordinate 
gives the value of D at which activation occurs at each site when it is reached by the activation 
fronts (Fi, continuous line and F2, dotted line). Each site is activated first by F\, and later by 
F2.The abscissa is the distance travelled by the fronts from the same arbitrary reference point 
(divided by L). Right panels: successive T act (x) (thick line) and T rep (x) (thin line) at each site. 
Each successive set of lines corresponds alternatively to the T ac t (x) and T rep (x) profiles associated 
to Fiand F2. The locations of F\ and F2 at any time t are the two positions where T ac t (x) = t. 
The successive values of D at one location are the differences between T ac t (x) and the preceding 



stimuli applied on period-1 SW reentry as a function of the timing of the pulses and the 
length of the loop. 
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A. Functional heterogeneity of refractoriness 

Complex dynamics can be induced by a second stimulus S2 thanks to the asymmetric 
profile in T rep (x) left by the interaction of first stimulus S\ with the reentry activation 
front F. When Apsi> the time between the last passage of F and the onset of Si, is 
beyond the vulnerable window, Si produces both a retrograde (Ri) and an antegrade (Ai) 
activation front. As illustrated in fig. |2K, the key factor determining the dynamics that 
can be induced by S2 is the region located between the site of stimulation (xs) and the 
site of the collision between F and R\ (x c , identified by the arrow in|2K). T rep (thin line 
in fig. [2JA.) is minimum near x$ and reaches its maximum at x c . The IL and ID models 
produce the same profile of T rep (x) (thin lines and triangles, respectively, in|2K), showing 
that the ID model that was initially developed to describe sustained reentry also provides an 
appropriate low-dimensional representation of the dynamics when stimulations are applied. 
The location of x c as well as the profile of T rep (x) around x$ depend on both Ap_si and L. 
Figure EJ3 shows the profile of T rep (x) obtained from loops of different lengths stimulated 
at the same diastolic interval Dsi(x s ) = A^_ S1 (L) — A F (L), where A F (L) is the action 
potential duration of the stable reentry for each L. The position of x c is shifted to the left 
(arrows in|23) because the collision is delayed on longer loops. However all the loops have 
the same invariant profile of T rep in the time and space interval that they share before the 
collision. If Si is applied at larger D S i value, the distance between xs and x c is shortened, 
A(xs) and T rep (xs) are increased, such that the extent and depth of the cusp in T rep around 
Xs are diminished. 

B. Initiating a second antegrade propagation 

The spatial profile of T rep for short A F _s\ is asymmetrical, with a sharp gradient between 
xs and xc, and a more gradual increase at the right of xs- Owing to this asymmetry, the 
outcome of S2 depends on Asi_s 2 , the time interval between the onset of the two stimuli. 
Figure H3 illustrates a case in which S2 is applied after the collision of F and R\, at an 
instant where Ai, the antegrade front created by Si, still has not reached x c . S 2 creates 
both an antegrade (A 2 ) and a retrograde (R 2 ) front, but R 2 is blocked between xs and x c . 



9 




B 



300 
200 
100 




<" - 
' ' \\ 

flu 

y \ 
y \ 






Z_=34cm 

- - L=30 cm 
L=25cm 



-10 -5 



5 

x—x (cm) 



10 



FIG. 2: A) T act (x)(thick line) and T rep (x) (thin line: IL model, >: ID model) following the 
application of a stimulus S\ at a site x s on a loop of 23 cm holding a period-1 SW reentry. S\ 
produces an antegrade front A\ propagating to the right, and a retrograde front R\ that collides 
with the front at the reentry at the site x c indicated by the arrow. B) T rep (x) around the stimulation 
site xs for loops of 25, 30 and 34 cm, using the ID model. Af-si was chosen such that D(xs) = 
— Ap(L) was the same on all loops, Ap(L) being the duration of the action potential of the 
reentry for each case. x c (|) is shifted to the left for longer L. The circle indicated the site where 

\dT re p/ dx\ Q 5 max- 
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FIG. 3: IL model, on a loop of 23 cm: S% creates both an antegrade (A2) and a retrograde front 
i?2-^2 is blocked between xs and x c shortly after its creation in the refractory tail left by the 
retrograde front i?iproduces by Si. 

Thereupon, the system is left with two antegrade fronts (Ai and A 2 ). This occurs as long 
as R2 does not propagate beyond x c to collide with A\, in which case A2 is left alone to 
perpetuate the reentry. This is an alternate scenario of unidirectional block that creates a 
propagating wave in the same direction as F and A\. 

All the complex dynamics occur in the range of Asi_s2 for which R2 is blocked between xs 
and x c whereas A 2 propagates. This depends on the profile of T rep (x) left by Si, which was 
shown to be an invariant function of D$i in fig. |2J3. Figure 0] shows the global characteristics 
of the dynamics in the [DsiA$i-s2 ] parameters plane for two values of L (L = 25cm < 2L min 
and 2L™°jf c_1 < L = 34 cm < 2L cr u). The parameter plane can be divided in three areas. 
In the region labeled at low A$i-S2 values, S2 is applied during the refractory period and 
does not produce a response. In the upper region, labeled " R 2 — ►", R2 propagates beyond 
x c , collides with Ai, and A2 is left alone to maintain the reentry. In the middle area, R2 
is blocked between x c and x s and complex dynamics may occur. The phase plane area in 
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which i?2 is blocked, for -Dsi, between and ~ 50 ms (dotted vertical line), is the same for 
the two L. The specific subsets in which complex dynamics occurs (represented by different 
shaded areas in fig. HJ) change with L, and are discussed later. 

The article is focused on the area where R2 is blocked, between D$i = and ~ 50 ms. The 
lower limit of the area coincides with A {Dsi) and is close to the action potential duration 
restitution curve A SP (D). The upper limit is nearly a constant, around Asi-s^max ~ 204 
ms. Appendix |X] shows that this A$i-S2 upper limit is set by the locus between x s and x c 
where \dT rep (x) /dx\ is equal to the maximum conduction time. This locus, indicated by 
a circle for the specific D$i illustrated in I2J3, does not depend on L, which explains why 
the upper limit for the block of R2 is identical in the two loops. \dT rep (x) jdx\ remains 
everywhere below the maximum conduction time if D$i is too long, explaining why the area 
of complex dynamics disappears beyond Dsi ~ 50 ms. It is also demonstrated in Appendix 
Klthat the maximum slope of the A SP (D) function must be greater than ~ 1 to allow the 

nnn 

block of i?2- The same condition that controls the stability of the period-1 reentry |20 U21LI27| 
thus determines if complex dynamics can be induced by 



C. Interactions between the two antegrade propagating fronts 

Once A2 has started to propagate and R2 has been blocked, there are 4 possible outcomes: 

1) Ax is blocked and A 2 perpetuates the SW reentry, 

2) A 2 is blocked, and Ai maintains the SW reentry, 

3) Ai and A 2 are blocked, and reentry is annihilated 

4) neither A\ nor A2 are blocked, and there is a transition to DW reentry. As seen in 
table HI this last option can only occur if L > 2L min . 

The next section discusses the cases 3) and 4) in which the system does not return to the 
original SW period-1 reentry. 



1. Reentry annihilation by alternans amplification (A\ and A2 are blocked) 

The three new scenarios of reentry annihilation reported in are shown in fig. These 
scenarios of termination differ in the number of revolution made by the A\ and /^activation 
fronts before they are blocked. Accordingly, we introduce the notation A™ H , meaning 
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FIG. 4: IL model: Outcome of the S\ — S2 protocol applied on a SW period-1 reentry in the 
[Dsi, A$i-S2] plane for two loops with L = 25 and 34 cm. The two loops have a common lower 
Region (0) where S2 is applied in the refractory period, and upper region ( R2 — >) where R2 
propagates beyond x c . In the area in between, R2 is blocked and complex dynamics may occur. 
They take place in a specific subset for each L: L = 25 cm, ligth gray, L = 34 cm, dark gray. 



that the front Aj = {1,2} is blocked (H) after m turns around the loop. Let consider the 
simplest case ([A\ H, A\ H], fig. in which both A\ and A 2 are blocked after one rotation. 
Ai is blocked first near xs, when it reaches the refractory tail left by A2 and i?2- This 
occurs because A\ has already completed a fraction of its rotation when S2 is applied. As 
a consequence, A\ comes back to reactivate the sites near xs after a time interval much 
shorter than the period of rotation T of the stable reentry. The block of A 2 takes place 
between x c and the site where R2 was blocked. When A\ travels in this zone, the region 
has last been excited by R\, and a time interval longer than T has elapsed since this last 
excitation. Besides, both the action potential and T rep associated to Ri were short since 
Si was premature. As a consequence, A\ produces action potentials that are longer than 
those of the stable reentry. Since the time between the passage of A\ and the return of 
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A 2 is also shorter than T, A 2 is blocked. The block of A 2 thus results from a process of 
amplified alternation in a region between xs and x c . The premature R\ that creates short 
action potential is followed by the late Ai generating long action potential. 

In the two other scenarios, A\ nor A 2 are blocked after either 2 or 3 turns (\A\ H, A\ H], 
\A\ ~\,A 2 H] in fig. [HJ3-C). Each passage of A\ leaves around xs a convex profile of T rep , 
which is turned to a concave profile by the subsequent passage of A 2 . The panel B and C 
illustrates the process of alternans amplification. Alternans amplification may end up by 
the blockade of one of the front, in which case propagation reverts to SW reentry. It may 
saturated, leading to DW reentry, or may lead to annihilation. The next section explores 
the conditions for annihilation. 

a. The scenarios of annihilation by alternans amplification occurs over limited range of 
L The conditions leading to each type of alternans amplification annihilation vary with 
L, and none exists for L >~ 30.5 cm. The four panels of fig. El picture the extent of the 
different zones of annihilation for L ranging from L cr u = 19.5 cm, the minimum L with 
stable period-1 SW reentry, to the limiting L w 30.5 cm value. At short L (L = 20.5 cm, fig. 
EK)j there is a large [D51, A5i_g 2 ] area with [A\ H, A\ H] block located at low D S \ values, 
and a small adjacent area of [A \ -\,A^ H] block. The zone of [Ai ~\,A 2 H] always remains 
minimal; and it is the first to disappear at L « 28.5 cm. The zone of [A\ H, A 2 H] appears 
at intermediate L, expands, and is the last to disappear. 

b. Disappearance of the [A\ ~\,A\ -\] annihilation Taken together, the maps of fig. 0] 
and fig. ElA-C shows that the area of [A\ H, A\ H] is embedded in the larger region, invariant 
with respect to L, in which R 2 is blocked between x c and xs- The zone of [A\ H, A\ H] is 
located at intermediate [\s\-S2-, just over the area, labeled [A\ — >,A\ H] (fig. HJ), in which 
Aican propagated beyond xs, but A 2 is stopped by the refractory tail left by A\. As L is 
increased, the lower boundary of [A\ H, A\ H] is lifted, thus diminishing the [A\ H, A\ H] 
area until it disappears completely. The lost of [A \ -\,Al, H] annihilation is thus caused by 
the inability of the system to block A\ at its first return. 

Appendix [B] provides the conditions needed for Ai to block in the tail of A 2 and proves 
that there is a limiting L beyond which this cannot happen. To summarize: 1) Increasing 
produces longer action potential for S 2 , delays T rep (xs), and thus augments the 
likelihood of A\ to be blocked. 2) However, the increase of Asi-S2 and T rep (xs) are bounded 
by the condition of R 2 being blocked between xs and x c , and these limiting values are 
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FIG. 5: IL model: profile of T ac t (thick line) and T rep (thin line) for the three different scenarios 
of annihilation caused alternans amplification model in the ring. A) L = 23 cm, D$i = 10 ms, 
= 160 ms, [A] -\,A\ H], B) L = 25 cm, D S i = 47 ms, A 5 i_ 52 = 190.4 ms, [A\ -\,A% H], 
C) L = 25 cm, Dsi = 47 ms, A5i_s2 = 185.1 ms [A\ H,^ ^i an d ^2 are pictured in panel A, 
but only 5*2 is shown in panel B and C. 

independent of L. 3) Since increasing L delays the return of Ai, there is a length L from 
which the return of A\ always occurs after the limiting T rep (xs) value. 

c. Opening the way to more than one rotation for both fronts Figure El shows that, for 
each value of L, the zone of [A\ H, A\ H] and \A 2 { 3 H, A% 3 H] blocks are disjoint, being 
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respectively located at low and high D$i- As seen in fig |Bp, the zones of [A 2 ^ H, A 2 ' 3 H] 
remain located at large Dsi even at values of L for which [A\ H, H] does not exists 
anymore. However, enlarging L also leads to the appearance and extension of areas in 
which A\ and A 2 persist together for an increasing number of turns before one of the front 
is blocked (up to 7 activations in HOD). As L is increased, the areas of the parameters plane 
associated to these other forms of transient complex dynamics extend toward low D$i, to 
finally cover all the range from D$i = to 50 ms when L = 37.5 cm, a value close to 2L crit 
from which DW reentry becomes possible. The areas of the parameters plane associated 
to these behaviors with coexistence of A\ and A 2 for multiple turns form a sequence of 
contiguous parallel tongues. 

All the higher modes of block and complex dynamics require that both A\ and A 2 travel 
beyond x s , at least at their first return. Conditions for the block of A 2 are derived in 
Appendix O based on an approximation using the A SP (D) equation and rules for A\ and 
A 2 return cycles. It shows that the block of A 2 depends on the balance between the return 
cycle of A\ and A 2 . 

2. Transition to double-wave reentry (neither A% nor A 2 is blocked) 

In both the ID and IL model, annihilation by alternans amplification is impossible from 
L > 30.5 cm. In fact, for 31 > L > 30.5 cm, only transient complex dynamics with a 
final return to period- 1 SW reentry are observed. However, the maximum number of turns 
with Ai and A 2 co-traveling grows, just as the number of tongues in the [Dsi, A-si-S2 ] 
plane associated to different numbers of turns during which the two fronts coexist. Each 
new tongue appears at low D S1 and A S1 _ S2 values, and expands as L is further increased. 
Finally, a first transition to QP mode-0 DW reentry is detected at L = 31 cm (fig. 
Transition to D W reentry thus appears as the asymptotic limit of the prolongation of the 
transient propagation with two fronts. However, transition to DW reentry begins much 
beyond L = 2L min = 25.6 cm, the value at which sustained mode-0 QP DW reentry starts 
to exist. In fact, at L — 31 cm, the system is rather in the range of L for which both DW 
mode-0 and mode-1 solutions coexists. 

For longer L, (fig. 03) the region with transition to mode-0 DW reentry expands, as 
it was the case for all the zones associated to two fronts transient propagation created at 
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shorter L. Transition to mode-1 DW reentry also appears. However, the transition to 
mode-1 appears at low Dsi, but at two disjoint intermediate Asi~s2 values, giving rise to 
two separated areas. 

For still longer L, the areas of mode-0 and mode-1 transition enlarge until covering, from 
L = 37.5 cm, all the area where R 2 is blocked. Hence, from this value of L, transition to 
DW reentry occurs in all instances where R 2 is blocked. Finally, beyond 2L cr i t = 39 cm, all 
transitions go to period-1 DW reentry, which is the only DW solution that remains. 

3. Selecting the mode of QP reentry 

Increasing the length of the pathway correlates with an increase in the complexity of 
the dynamics and to the transition to sustained DW reentry. The results presented in fig. 
13 show that changing the timing of the stimuli when 2L min < L < 2L crit can induce the 
transition either to the mode-0 or mode-1 DW QP reentry. The IL model was simulated with 
L = 37.5 cm (2L m [ n < L < 1L cr it but near 2L cr n) to circumscribe the basin of attraction 
in the parameter space associated to each DW QP reentry. The large area in which two 
antegrade fronts are created (Dsi < 50 ms) is separated between two regions, with lower 
Asi-S2 converging to mode-0, and higher A$i-S2 to mode-1. In this last section, we compare 
the transient dynamics leading to each of two modes of DW QP reentry. 

Figure |H1 shows the spatial profile of D associated to A\ and A 2 for two cases converging 
respectively to mode-0 (panel A) and mode-1 (panel B). In this representation, the passage 
of A 2 (Da2 { z )i dotted line) at each location is followed by that of A x (Dai (z) ,full line). 
At first, A2 propagates with a short D, producing brief action potentials everywhere along 
the loop. As a consequence, A\ travelling afterwards meets long D (z). Dai (z) has its first 
maximum (Pi) when it travels between x c and the point where R 2 was blocked, such that 
Da2 has a minimum at the same location from which it increases slowly until its next return 
in the same region. From there, Dai (z) starts to alternate between long and short values 
with a spatial period close to 2 turns, and D^i (z) follows a complementary profile. This first 
phase of the propagation, lasting for approximately 8 turns, can be labeled as concordant 
alternans since Da2 and Dai have alternating values and that each remains short or long 
for at least one complete rotation. This pattern does not correspond neither to mode-0 nor 
mode-1, since both stabilized DW QP solutions have a wavelength less than L. 
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However, from the beginning (k = 1,2), Dai already shows a second spatial oscillation 
in D that is superimposed to the concordant alternans. The structure of this oscillation, 
which has a wavelength close to one turn, makes the difference between the cases converging 
to mode-0 and mode-1. In the former case, the oscillation embeds two peaks {Pi, P2} (-Dai 
for mode-0, continuous line in fig. UK), while in the later case, it has three peaks {Pi, 
P 2 , P 3 } (Dai for mode-1, continuous line in fig. |HJ3). These superimposed spatial variations 
persist while the concordant alternans dissipate. During this process, the position of the 
peaks does not change much. As the amplitude of the concordant alternans decreases, the 
respective height of the discordant alternans increases up to a point where the boundary 
with large gradient in D begins to move around the loop due to the quasiperiodic nature of 
the propagation. This travelling mechanism is akin to the propagation of paced discordant 



alternans on a cable of cardiac tissue 
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The main difference between the transition to mode-0 and mode-1 is the presence of 
the third peak P 3 in Dai- The P3 peak of Dai is induced by the increase of Asi-S2- On 
one hand, a larger Asi_s 2 produces a longer T rep such that A\ meets a lower D and the 
amplitude of Pi is reduced. But a larger A5i_52 also implies that A2 travels faster, comes 
back sooner to the stimulation sites, and set the stage for a new maximum in the Dai profile. 

Of course the transition from two to three peaks is continuous process since a similar 
variation of D with less amplitude already exists with the transition to mode-0 in fig. |HK. 
It means that there must be a minimal spatial profile that corresponds to the boundary 
between the two basins of attraction (the transition to either mode-0 or mode-1). 

IV. DISCUSSION 

Alternans amplification, leading either to reentry annihilation or transition to DW reen- 
try, can exist if the slope of the APD restitution curve is larger than one. Hence, the 
condition on the slope of the APD restitution curve that is mandatory for the existence of 
sustained QP reentry|2l| also enables double-pulse stimulation to produce a new mode of 
unidirectional block, in which only an antegrade front propagates away from the stimula- 
tion site. On any closed circuit with two activation fronts travelling in opposite direction, 
this opens a large spatio-temporal window in which an ectopic focus or an external source 
firing twice can start a reentry. This is consistent with the use of burst pacing as a standard 



18 



experimental and clinical procedure to start tachycardia 2^ 13 |30|> |3l| . In a loop already 
holding a SW reentry, the Si — S 2 protocol can induced the simultaneous propagation of two 
antegrade fronts whose final outcome depends on the timing of the stimuli and the length 
of the loop. The transient or persistent coexistence of two antegrade fronts is a new type 
of dynamics in which the effects of the stimuli cannot be represent as perturbations of a 
imit-cycle, as it has been done for models in which the steep slope criteria was not fulfill 



aalaai- 

In clinical and experimental investigations, stimulations are currently used to either st ud) 



utner st udy 
lias 0, Q|. 



the characteristics of reentry circuits through resetting or to stop the tachycardi 
Most often, stimulations are applied at one site and propagation is assessed through one 
or a few recording electrodes. We may consider the dynamics that will be observed with 
this setting for each of the three scenarios of annihilation by alternans amplification. Since 
the distance traveled by -Ri from x s to x c usually covers a small portion of the reentry 
pathway, electrodes are much more likely to be positioned in the complementary portion of 
the circuit. For [A \ ~\,A\ H], such electrodes will detect the last passage of F before the 
collision with Rx, then Ax and A 2 . The detection of A\ and A2 would clearly exclude classical 
unidirectional block being responsible for the annihilation. In fact, this modeling study was 
initiated after a set of experimental and clinical studies on flutter using multichannel (4 to 8 
channels) recordings 3]. In these, cases of annihilation were reported in which F, A\ , A 2 
were detected, and in which the propagation of A\ and A2 was blocked in the segment of the 
circuit were the collision of F and A% was presumed to have occurred. This scenario, that 
was called collision block, is consistent with the [A \ ~\,A\ H] block by alternans amplification 
[3]. In the case of [A\ ~\,A 2 H] block, A\ and A 2 would be detected twice (see 03). Since 
the time intervals from each A\ to the next A 2 are longer than that between A 2 and the 
following Ax, the sequence of the differences between activation times alternates around a 
value shorter than the period of the original reentry, the block occurring after the longer 
interval of the time series. For [A\ H, A\ H], each front is seen thrice, with similar oscillation 
in the time series of the difference. However, the structure of oscillation of the cycle as well 
as the last value before annihilation depend on the position of the electrode in the circuit. In 
a protocol for annihilation by unidirectional block, the detection of A\ and A 2 would indicate 
that the stimuli were beyond the vulnerable window and would trigger the application of 
a new stimulus at shorter coupling interval. However, this could rather reinitiate a reentry 
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bound to stop by alternans amplification. 

To allow the unidirectional block of R2, Si must be applied early beyond the vulnerable 
window, in t he p ortion of the excitable gap referred as partially refractory by electrophys- 
iologists [ij The prematurity of Si, in conjunction with the high slope of the APD 
restitution curve and of the 9 dispersion curve, create a concave asymmetrical profile of 
T rep (x) around x s . For a given Si, the Asi-S2 intervals for which R 2 is blocked depends 
on the local dynamics around x s that is again determined by the APD restitution and the 
dispersion of 9. In fact, when the parameters D$i is used to describe the timing of the first 
stimulus, the [D$i, Asi~s2] range to get a block of R2 becomes independent of L. 

The pivotal role of the A SP (D) and 9(D) functions is further confirmed by the capacity 
of the ID model to reproduce the dynamics of the ionic model. However, in order to avoid 
discontinuity in T rep (x) in the region where R2 is blocked and get a correct representation of 
the dynamics, spatial averaging must be included in the computation of the spatial profile of 
T rep (x). Originally, spatial averaging was added to the ID model to reproduce the modulation 
of the repolarization by resistive coupling in order to correct the shortcomings of the model 
regarding the details of the bifurcation from periodic to QP propagation [t| . Including spatial 
averaging becomes even more essential when stimulations and blocks produce steep gradients 
in T rep (x) Q- 

While the wide range of [D S i, A S i-s2] for which R 2 is blocked does not depends on L, 
this is not the case for the subset of this interval leading to reentry annihilation by alternans 
amplification. As shown in fig. the interval of block is very large for L ~ L cr i t and 
decreases gradually until it disappears at L ~ 2L™?f e " 1 . For the MBR model used herein, 
annihilation by alternans amplification exists for 30.5 > L > 19.4 cm, and it occurs on a 
range of [Dsi, Asi-S2\ that remains much wider than the ~ 1 ms standard vulnerable window 
for most of this range of L. This is consistent with the result of a clinical study in which dual 
pulses stimulation was found to be four times more effective than single stimulus to stop 
monomorphic ventricular tachycardia in man 35]. As L is increased, the [Dsi, Asis2] area 
of R2 block also encloses a growing number of tongues with increasingly prolong coexistence 
of Ai and A 2 , a process that culminates in the appearance of zones of transition to DW 
reentry. Then the zones with transition to DW reentry extend until they cover completely the 
area in which R 2 is blocked. From L ~ 2L crit , blocking R 2 yields automatically to sustained 
DW reentry. Hence, transition to DW reentry does not necessarily imply a heterogeneous 
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substrate as it is been proposed elsewhere 



361 l37l | , but can also be achieved in a homogeneous 



medium through the creation of a functional heterogeneity by a limited number of electrical 
stimulations. Annihilation and transition to DW occur on separated ranges of loop lengths. 
On shorter loop, reentry annihilation is produced when alternans amplification reaches an 
amplitude high enough to block y^and A 2 . On longer loop, the distributed alternans saturate 
at an amplitude that still permits sustained propagation. In our version of the BR model, the 
transition from SW QP to period-1 reentry at L CT u is supercritical. However, we have shown 
that, for other sets of parameters, the bifurcation is subcritical, with bistability between QP 
and period-1 reentry near L cr n In these cases, it is possible that the two-stimulations 
protocol applied to the SW period-1 reentry near L crit could induced a transition to SW 
QP propagation, a phenomenon that was not possible with the instance of the MBR model 
used in this paper. 

Our results are consistent with different clinical and experimental observations, and open 
the possibility to design more effective anti-arrhythmic pacing strategies. However, the 
modeling studies must be extended to more realistic representations to evaluate properly 
potential applications. Preliminary results from an ongoing work on a two-dimensional 
annulus show that block by alternans amplification can still be obtained on this setting, but 
that additional scenarios are possible, including termination through transient fibrillation 
that has also been observed in real cardiac tissue. Tissue heterogeneity, either at the level 
of the ionic properties or of the cells coupling, could also be important since termination of 
reentry has been obtained in a one- dimensional loop model embedding a small area of slower 
conduction but using an ionic model with minimal APD restitution properties 

HQ. it 

remains to be seen if annihilation based on alternans amplification would be amplified or 
reduced by the inclusion of spatial inhomogeneity. Investigation will also have to be extended 
to bidomain model in order to get a more proper representation of the stimulus. The 
simplified representation used herein can be an acceptable approximation for low amplitude 
stimuli or to mimic the effect spontaneous firing of a group of cells. However it is known 
that current spread of the stimulus depend on the properties of the external and internal 
medium [iol . \41 \. The MBR representation of the ionic properties is also oversimplified. 
However, since most phenomena described in this paper occurs in the few first beat after 
the stimulations and can be explained from the APD restitution and speed dispersion, it is 
unlikely that slow memory effects appearing in the dynamics of more complex model would 
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change the behaviour. Characterizing the APD restitution and speed dispersion of the more 
complex model in the range of frequency of repetitive activity associated to reentry should 
allow a prediction of the possible dynamics. 

V. CONCLUSION 

This work is a further illustration of the richness and diversity of the dynamics that 
can results from the restitution of APD and dispersion of the speed even in a simplified 
model of the tissue. It has revealed some unexpected behaviours, like block by alternans 
amplification, which can be much more prevalent than the mechanism of unidirectional block 
that are usually assumed to be dominant. The low-dimensional model whose behaviour is 
equivalent to the ionic model, provides a generic understanding of the dynamics that can 
used a guideline to investigate the effects of future complexification of the model. 

Acknowledgments 

This work was supported by grants from the Natural Sciences and Engineering Research 
Council of Canada (AV), the Fonds Quebecois de la Recherche sur la Nature et les Tech- 
nologies (PC), as well as by the technical and computer resources of the Reseau Quebecois 
de Calcul de Haute Performance. 



[1] G. R. Mines, Trans Roy Soc Can 4, 43 (1914). 

[2] L. H. Frame and M. B. Simson, Circulation 78, 1277 (1988). 

[3] L. H. Frame and E. K. Rhee, Circ Res 68, 493 (1991). 

[4] J. M. Pinto, J.N. Graziano, and P. A. Boyden, J Cardiovasc Electrophysiol 4, 672 (1993). 

[5] E. Jalil, B. Mensour, A. Vinet, and T. Kus, Can J Cardiol 19, 244 (2003). 

[6] M. Courtemanche, L. Glass, and J. P. Keener, Phys Rev Lett 70, 2182 (1993). 

[7] A. Vinet and F. A. Roberge, Ann Biomed Eng 22, 568 (1994). 

[8] A. Vinet, Journal of Biological Systems 7, 451 (1999). 

[9] A. Vinet, Ann Biomed Eng 28, 704 (2000). 



22 



[10] A. Garfinkel, Y. H. Kim, O. Voroshilovsky, Z. Qu, J. R. Kil, M. H. Lee, H. S. 

Karagueuzian, J. N. Weiss, and P. S. Chen, Proc Natl Acad Sci 97, 6061 (2000). 
[11] Z. Qu, A. Garfinkel, P. S. Chen, and J. N. Weiss, Circulation 102, 1664 (2000). 
[12] P. Delia Bella, G. Marenzi, C. Tondo, D. Cardinale, F. Giraldi, G. Lauri, and 

M. Guazzi, Am J Cardiol 68, 492 (1991). 
[13] M. Heldal and O. M. Orning, Eur Heart J 14, 421 (1993). 

[14] A. Prakash, S. Saksena, M. Hill, R. B. Krol, A. N. Munsif, I. Giorgberidze, P. Mathew, 

and R. Mehra, J Am Coll Cardiol 29, 1007 (1997). 
[15] P. Comtois and A. Vinet, Chaos 12, 903 (2002). 

[16] W. L. Quan and Y. Rudy, Pacing Clin Electrophysiol 14, 1700 (1991). 
[17] R. M. Shaw and Y. Rudy, J Cardiovasc Electrophysiol 6, 115 (1995). 
[18] H. Fei, M. S. Hanna, and L. H. Frame, Circulation 94, 2268 (1996). 
[19] A. Vinet and F. A. Roberge, J Theor Biol 170, 183 (1994b). 

[20] M. Courtemanche, J. P. Keener, and L. Glass, Siam Journal on Applied Mathematics 
56, 119 (1996). 

[21] P. Comtois and A. Vinet, accepted for publication in Phys. Rev. E (2003) 

[22] M. A. Watanabe, F. H. Fenton, S. J. Evans, H. M. Hastings, and A. Karma, J 

Cardiovasc Electrophysiol 12, 196 (2001). 
[23] J. J. Fox, J. L. McHarg, and R. F. . J. Gilmour, Am J Physiol Heart Circ Physiol 282, 

H516 (2002). 

[24] B. Echebarria and A. Karma, Phys Rev Lett 88, 208101 (2002). 

[25] J. Starobin, Y. I. Zilberter, and C. F. Starmer, Physica D 70, 321 (1994). 

[26] B. Mensour, E. Jalil, A. Vinet, and T. Kus, Pacing Clin Electrophysiol 23, 1200 (2000). 

[27] E. Cytrynbaum and J. P. Keener, Chaos 12, 788 (2002). 

[28] A. Vinet, R. Cardinal, P. LeFranc, F. Helie, P. Rocque, T. Kus, and P. Page, 

Circulation 93, 1845 (1996). 
[29] J. M. Cao, Z. Qu, Y. H. Kim, T. J. Wu, A. Garfinkel, J. N. Weiss, H. S. Karagueuzian, 

and P. S. Chen, Circ Res 84, 1318 (1999). 
[30] F. Helie, A. Vinet, and R. Cardinal, J Cardiovasc Electrophysiol 11, 531 (2000). 
[31] S. Nattel, D. Li, and L. Yue, Annu Rev Physiol 62, 51 (2000). 
[32] T. Nomura and L. Glass, Physical Review E 53, 6353 (1996). 

23 



[33] L. Glass, Y. Nagai, K. Hall, M. Talajic, and S. Nattel, Phys Rev E 65, 021908 (2002). 

[34] A. Heisel, J. Jung, M. Stopp, and H. Schieffer, Eur Heart J 18, 866 (1997). 

[35] J. M. Almendral, M. E. Rosenthal, N. J. Stamato, F. E. Marchlinski, A. E. Buxton, 

L. H. Frame, J. M. Miller, and M. E. Josephson, J Am Coll Cardiol 8, 294 (1986). 
[36] J. Brugada, P. Brugada, L. Boersma, L. Mont, C. Kirchhof, H. J. Wellens, and 

M. A. Allessie, Circulation 83, 1621 (1991). 
[37] J. Cheng and M. M. Scheinman, Circulation 97, 1589 (1998). 
[38] S. Sinha, K. M. Stein, and D. J. Christini, Chaos 12, 893 (2002). 
[39] S. Sinha and D. J. Christini, Phys Rev E 66, 061903 (2002). 

[40] A. E. Lindblom, B. J. Roth, and N. A. Trayanova, J Cardiovasc Electrophysiol 11, 274 
(2000). 

[41] J. P. Keener and E. Cytrynbaum, J Theor Biol 223, 233 (2003). 

APPENDIX A: THE UPPER LIMIT OF THE A 5 i_ 52 INTERVAL FOR DOUBLE- 
WAVE CREATION 

Three conditions must be fulfilled for the block to occur: 1) R 2 , the retrograde front 
created by S 2 , must be blocked between the stimulation site x s and x c , the locus of the 
collision between the reentry front F and retrograde front R± created by Si; 2) afterward, 
Ai, the antegrade front produced by Si, must be blocked when it returns near x s ; 3) finally, 
A 2 , the antegrade front produced by S 2 , must also be blocked when it travels between x s and 
x c The following three appendixes formulate the constraints associated with each of these 
conditions. In all three appendixes, we consider that F, ^and A 2 travels toward increasing 
value of x, we use T rePjF (x s ), the repolarization time associated to the passage of the last 
activation front of the reentry F at xs, as the reference time t — 0, and introduce a spatial 
coordinate y = (x s — x) > to follow the retrograde fronts Ri and R 2 . 

R 2 can propagate as long as its activation time T acty R 2 (x) is larger than T rePj m(x), the 
repolarization time associated to Ri. In the limiting case T ac i fi 2 {oc) — T re p R i{x) j Dr 2 {x} — 
T a ct,R2 (^O T rep ,ri(x) = 0, and R 2 travels with the maximum conduction time -dT actjR2 /dx = 
c t (0) =c tjmax . When R 2 reaches a point x c < x < x s with D R2 (x ) = 0, it can continue to 
propagate if \dT rePj Ri(xo) / dx\ < Q imax . The limiting case for the propagation of R 2 is thus 
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T re p,Rl(Xo) = T act ^ R2 (x ) (Al) 



dx 

We neglect the effect of coupling in the calculation of T rePj m to obtain 

T rep ,m (y) ^ T actiR1 (y) + A SP (D R1 (y)) (A2) 

where Dm (y) — T act> m(y) — T rep ^(y) is the diastolic interval associated to the propagation 
of R\. From our choice of reference time and the definition of y, T actt m(ty — Asi- I n the 
MBR model, the c t dispersion curve is very steep such that the conduction time is minimal, 
except for a short interval of D close to 0. As a consequence, we approximate that both 
Ri and F have been propagating with the minimum conduction time (i.e. maximum speed) 
Q,min, such that 

T re p J F(y S ) — Ct,min?/ 

D m(y) = T act)R1 (y) ~ T reP)F (y) ~ D S1 + 2c tMn y 



Substituting these relations in eq. IA2I yields 

(y) ~ D S i + Ct Mn y + A SP (D S i + 2c tMa y) (A3) 
Taking the spatial derivatives of eq. IA3I yields 

dT rep> m I - - dA SP 



Ct,min I 1 "I - 2 



dy \ dD 

which, thanks to eq. IA3[ gives 



D m (y), 



dA SP 



Ct ,max 

D m {yo)=D^ 2c t;min 



dD 

The existence and location of the critical point depends on the slope of the restitution curve, 
and on the relative difference between the maximum and minimum conduction time. From 
the c t (D) given by eq. H, Q imax ~ 3.13Q imm , and 

dA SP 



dD 



~ 1.065. (A4) 

ffll(9o) = max 
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Using eq. |3] to solve this equation, we obtain that D mSuX ps Dsi + 2ct )IO myo — 70 ms. This 
result means that, if Dgi > D max ~ 70 ms, \dT rePt m/dx\ < c i)max everywhere between x c 
and x s , and R2 cannot be blocked. It means also that if D$i < D max , there is a critical 
point x c < x < x s whose position depends on D$i but is independent of L. The last step is 
to determine the maximum A si _s 2 (A si _s 2 ,max) to get a block of R 2 at the x . Assuming 
again that i? 2 travels with the minimum conduction time, the first condition of eq. IA1I 
becomes 

D S i + c timin y + A (D max ) = D S1 + Asi 

—52, max Q,minZ/ 

reducing to the limit 

A ASP ( p) \ 

Si— S2' m ax V max/ 

The maximum A5i_s 2 for the block of R2 is independent of L and of D$i, provided that 
■ Dsi < -Dmax- For S2 to induce propagation Asx.^, must also be > APD(Dsi), which is 
the refractory period at x s . In summary, for all < Dg\ < -D m ax — 70 ms, R2 blocks if 

APD SP (D S1 ) < A S1 _ S2 < A SP (-D m ax) ■ 

The main approximation used herein is that both Ri and R2 propagate everywhere with 
the minimal propagation time. With regard to R±, the error introduces by the approximation 
is minimal unless D$i — 0. For R2, the conduction time is obviously underestimated when 
its activation time comes close to T rePt m. The approximation thus overestimates -D max , whose 
value is around 50 ms for the ID and IL model, compared to the -D max = 70 ms provided 
by the approximation. 



APPENDIX B: BLOCK OF A 1 AFTER ONE ROTATION ON THE LOOP 

Appendix |A] shows that VD51 G [0, -D max ], R2 exist and is blocked if A S1S2 e 
[APD SP (Dsi), A SP (Dmax)], a set of conditions that do not depend on L. The next event 
is the annihilation A\ when it returns near x s and hits the refractory tail left A2. Using 
T rep ^(x s ) as a reference time A\ comes back to x s at the time D$\ +T41, in which is the 
time taken by A\ to propagates over the loop on its first turn. It is blocked if the system is 
still refractory, which means 

T re p,S2 > D S 1 + T Al (Bl) 
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Neglecting the effect of coupling on repolarization, T reP) s2 is approximated as: 

T rep ,s2 ^ Dsi + A 51 _ 52 + A SP {D S2 ) (B2) 
in which -D52, the diastolic interval associated to 5*2, is estimated by 

D S2 ~ A 51 _ S2 - A SP (D S1 ) (B3) 
Substituting eq. IB2I and lB3l in eq. IB1I yields 

A S1 _ 52 + A SP (A S1 ^ S2 - A SP (D S1 )) > TAl , (B4) 

in which r Al is the only non- local term. For most ionic models and experimental prepara- 
tions, tax is a monotonic decreasing function of D$i (as in fig. 4 of ref. ^J). The effect 
of the prematurity on the return cycle comes from the steepness of the c t dispersion curve 
(eq. EJ) which is close to c tim i n as soon as D is greater than a few tenths of ms. Therefore, 
the prolongation of the return cycle depends on the limited region beyond x s over which A\ 
does not propagate at the maximum speed. Hence, we write 

tai (D S1 , L) w L Q , min + / (D S i) (B5) 

where / is maximum for D$i = and 

df 

< 0, 



dD sl " ' 
lim / = 

D S1 — >oo 

The condition IB4I for the block of A\ becomes 

A 5 i- S2 + A SP (A S1 _ S2 - A SP (D S1 )) > Lc tMn + f (D S1 ) . (B6) 

Consider 

A^ S2 + A SP (A^ S2 - A SP (D S1 )) = Lc t , min + f (D S1 ) (B7) 

as the minimum value of Asi_s2 to get a block of A\. 

Lets consider next the case where Dsi < D max and Asi_s2 = A SP (Dsi), for which 
Ds2 = 0. The condition for the block of Ai becomes 

A SP (D S1 ) + A SP (0) = Lc tMn + f (D S1 ) 
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which cannot be fulfilled since both D$i and A SP (0) are smaller than the diastolic interval 
of the free reentry. The block of A\ will occur from the minimum I\s\-S2 > A SP (Dsi)- 
This explains why the dynamical regime [A\ — A\ H] is found in the lower portion of the 
[Dsi, A S1 _52] area in which R 2 blocked, as seen in figure |U 

The condition for the block of R 2 are that Dgi G [0, -D max ] and A gl _ S2 G 
[APD SP (Dsi), A SP (D max )]. Since the right hand side of eq. IB 71 is a growing function 
of L, A^™_ S2 for any fixed value of Dgi must also increase with L until reaching A (D max ). 
Hence, there is a limit value of L from which A\ cannot be blocked on its first return. For 
Dg™(L), the value of D S i for which A l ^_ S2 = APD SP (D max ) = A 51 _ 52jmax , eq. IB7I becomes 

A 5 l-52, m ax + ^ P (A 51 _ 5 2,max - A SP (D l ™(L)) = L Q , min + / (D%?{L)) 

Shortening D l ™(L) increases both the left and right side of the equation, such that the 
value of Dg±(L) depends on the balance between the slope of df /dD (i.e. the change in the 
return cycle) and dA SP /dD (the change in repolarisation time at x s ). The derivative of the 
equation with respect to L gives: 

/ df , , \ dD'™ _ 

V dD sl aD si a n S2 ) dL " c *--n 

where a' Do = dA SP /dD\ D . In the MBR model, \df/dD sl \ < 1, whereas dA SP /dD > 1 
at low D value, such that dD x ^/dL < 0, as observed in the numerical simulations. 

APPENDIX C: BLOCK OF A 2 AFTER ONE ROTATION ON THE LOOP 

Appendix IB"1 shows that A\ is blocked only over a subset of the [Dgi, A S1 _ S2 ] area for 
which R 2 exist and is blocked (i.e. D S1 e [0,£W], A S1 _ S2 G [APD SP {D SI ), A SP {D max )}). 
With respect to A S1 _ S2 , the lower bound of this subset is greater than APD SP (D S i) and 
increases with L, while the upper bound remains constant at A SP (D max ). With respect to 
D s i,the variation of the limits of this subset is given by a complex expression that depends on 
the slope of both the restitution and dispersion curves. For the MBR model, the upper D$i 
of the subset decreases toward as L is increased. The last step is to obtain the conditions 
for the block of A 2 . 

A 2 blocks between x c and the locus where R 2 has stopped when it hits the refractory tail 
left by A%. Using the y coordinate, the condition for the block of A 2 is that there is a point 
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where 



T re p,Al (V) > T act,A2 (jj) (CI) 

T rep ,Ai(y) is approximated as T actjA1 (y) + A SP (D Al (y) , where D A1 (y) = T actA1 (y) - 
T re p,m (y) since R 2 has not propagated in this region. Using eq. IA3I of appendix |X] for 
Tre P ,m{y) and eq. |B5]of appendix |B] for T act)A1 (y) yields 

T repo i, A1 = D S1 + Lc tMn + / (D S i) - CtMnV + A SP {D Al (y) (C2) 
D A1 (y) = Lc tMn + / (D S1 ) - 2c tMn y - A SP (D S1 + 2c tMn y) 

T ac t,A2 (y) = Dsi + A-si-si + T A2(y) where r A2 is the time needed for A 2 to travel from the 
stimulation site to y. For A±, we have assumed in El that the prolongation of the return 
cycle was occurring mainly in a short region around x s beyond which A\ was travelling 
at maximal speed. The situation is different for A 2 . As it could be seen in fig. El A, the 
D associated to A 2 as it propagates away form x s , meaning that it speed of propagation 
diminishes. Nevertheless, we write 

r A 2(y) = Lc t , m - m + g (D S iD s2 ) - c t , min y (C3) 

With eq. ElandEHl the condition O for the block A 2 after one rotation becomes 

/ (D S1 ) + A SP (D A1 (y)) > A sl s2 + 9 (D S i, D S2 ) (C4) 

To get the [A\ H, A\ H] block, the conditions given by eq. IB6I and eq. IC4I must both be 
fulfilled, leading to the final condition: 

A SP {A S1 _ S2 - A SP (D S1 )) > Lc tMn + g (D S1 , D S2 ) - A SP (D A1 (y)) . 

with the supplementary constraint that A$i-S2 and D$i must remain in the interval where 
R 2 is block. Since both Asi_s 2 and Dsi are bound, there is a limiting L value beyond 
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which [A\ H, A\ H] block cannot occur. For the IL and ID model, we found that A2 was 
always stop when A\ was blocked, meaning that eq. IC4I was fulfilled whenever eq. IB6I 
was satisfied. However, the condition depends on the restitution and dispersion curve and, 
through g (D$i, -D52), the interaction of A 2 with the spatial profile of T repo i left by A\. 
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FIG. 6: IL model: Area with complex dynamics and alternans amplification in the (Dsi, A,si_s2) 
plane with A to D) L = 20.5, 25.0, 27.5, 28.5 cm. The gray scale indicates the number of turns 
during which A\ and Ai persist together. The regions with different annihilation are specified 
accordingly. Otherwise, one front ends up being blocked and propagation reverts to period- 1 SW 
reentry. 
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FIG. 7: IL model: area of complex dynamics in the D$i vs. Asi_s2 plane for L = 31 cm (panel 
A) and L = 34 cm (panel B). The gray code indicates the number of turns during which A\ and 
A-2 persist together. Arrows pinpoint the zones with transition to mode-0 and mode-1 DW QP 
dynamics. 
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